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Jet quenching theory using perturbative QCD is extended to include (1) elastic as well as (2) 
inelastic parton energy losses and (3) jet path length fluctuations. The extended theory is applied to 
non-photonic single electron production in central Au-|-Au collisions at = 200 AGeV. The three 
effects combine to significantly reduce the discrepancy between theory and current data without 
violating the global entropy bounds from multiplicity and elliptic flow data. We also check for 
consistency with the pion suppression data out to 20 GeV. Fluctuations of the jet path lengths 
in realisitic geometry and the difference between the widths of flnctuations of elastic and inelastic 
energy loss are essential to take into account. 


PACS numbers: 12.38.Mh; 24.85.-|-p; 25.75.-q 

Light quark and gluon jet quenching observed via 
TT, rj suppression [l| in Cu-|-Cu and Au-|-Au collisions at 
\/s = 62 — 200 AGeV at the Relativistic Heavy Ion Gol- 
lider (RHIC) has been remarkably consistent thus far 
with predictions 0-0. However, recent non-photonic 
single electron data B a i i 0 (which present an 
indirect probe of heavy quark energy loss) have signifi¬ 
cantly challenged the underlying assumptions of the jet 
tomography theory (see [HI). A much larger suppression 
of electrons than predicted was observed in the pr ~ 4—8 
GeV region (see Fig. [T]). These data falsify the assump¬ 
tion that heavy quark quenching is dominated by radia¬ 
tive energy loss when the bulk QCD matter parton den¬ 
sity is constrained by the observed dN/dy « 1000 rapid¬ 
ity density of produced hadrons. 

The observed “perfect fluidity” [13, of the sQGP at 
long wavelengths (j>t < 2 GeV) provides direct evidence 
for highly nonperturbative bulk dynamics [ij, [13] • Due 
to asymptotic freedom, a breakdown of perfect fluidity 
and nonperturbative effects are expected at pp several 
times greater than the mean thermal energy, 3T ~ 1 — 2 
GeV. Prior to these electron data, pQCD based jet 
quenching theory provided increasi ngly reliable predic¬ 
tions above pr > 5 — 7 GeV for the nuclear 

modification of light parton jets [3|-0- However, the 
non-photonic single electron data however raise the ques¬ 
tion of whether the novel nonperturbative physics of the 
strongly interacting Quark Gluon Plasma (sQGP) [3 
produced at RHIC could persist down to much smaller 
wavelengths than previously expected. This question is 
also of pragmatic importance because high pp jets can 
be utilized as calibrated “external” tomographic probes 
of the bulk sQGP matter only if their dynamics can be 
predicted reliably. 

The upper band of Fig. [T] shows that the predictions 
from considerably underestimate the electron nu¬ 
clear modification of data even out to pr ~ 8 GeV. This 
discrepancy points to one or more of (1) missing per¬ 
turbative QCD physics, (2) incomplete understanding of 
the initial heavy quark production and/or (3) novel non- 



FIG. 1: The suppression factor, Raa{pt), of non¬ 
photonic electrons from decay of quenched heavy quark 
(c-|-b) jets is compared to PHENIX 0 and STAR 0 
data in central Au-|-Au reactions at 200 AGeV. Shaded 
bars indicate systematic errors, while thin error lines 
indicate statistical ones. All calculations assume initial 
dNg/dy = 1000. The upper yellow band from m takes 
into account radiative energy loss only, using a fixed 
L = 6 fm; the lower yellow band is our new predic¬ 
tion, including both elastic and inelastic energy losses 
as well as jet path length fluctuations. The bands pro¬ 
vide a rough estimate of uncertainties from the leading 
log approximation for elastic energy loss. The dashed 
curves illustrate the lower extreme of the uncertainty 
from production, by showing the electron suppression 
after both inelastic and elastic energy loss with bottom 
quark jets neglected. 


perturbative mechanisms affecting partonic physics out 
to pt > 10 GeV. We note that pr ~ 8 GeV (single 
non-photonic) electrons originate in our calculations from 
the fragmentation and decay of both charm and bottom 
quarks with transverse momenta px ~ 12 ± 4 GeV (see 
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Fig. [5]in [HI)- 

Possibility (3) is of course the most radical and would 
imply the persistence of non-perturbative physics in 
the sQGP down to extremely short wavelengths. Pro¬ 
cesses can be postulated to improve the fit to the 
data but at the price of losing theoretical control 
of the to mog raphic information from jet quenching data. 
DGVW showed that by arbitrarily increasing the 
initial sQGP densities to unphysical dNg/dy^ 4000, the 
non-photonic electrons from heavy quarks can be arti¬ 
ficially suppressed to Raa ~ 0.5 ± 0.1. Thus, to ap¬ 
proach the electron data, conventional radiative energy 
loss requires either a violation of bulk entropy bounds 
or nonperturbatively large extrapolations of the the¬ 
ory. Even by ignoring the bottom contribution. Ref. IlSlI 
found that a similarly excessive transport coefficient [20|| , 
Qeff ~ 14 GeV^/fm, was necessary to approach the level 
of suppression of electrons in the data. 

Bottom quark jets are very weakly quenched by ra¬ 
diative energy loss. Using the FONLL production cross- 
sections, their contribution significantly reduces the sin¬ 
gle electron suppression m compared to that of the 
charm jets alone. The ratio Raa is not sensitive to the 
scaling of all cross-sections by a constant. However, it is 
sensitive to any uncertainty in the relative contribution of 
charm and bottom jets to the electrons [l^. Recent data 
from STAR on electrons from p-l-p collisions Q may in¬ 
dicate an even larger uncertainty in the production than 
expected from FONLL. However, PHENIX p-l-p to elec¬ 
tron data are compatible with the upper limit of FONLL 
predictions [2lL l22l |. similar to the comparison between 
FONLL and Tevatron data. 

The discrepancy between the ‘DGLV Rad only’ predic¬ 
tions and the data in Fig. [T]and recent work [3, 
motivated us to revisit the assumption that pQCD elas¬ 
tic energy loss [2^ is negligible compared to radiative. In 
earlier studies, the elastic energy loss [2^, [13, [Hi HI, HO, 
[3l| was found to be dE^^/dx ~ 0.3 — 0.5 GeV/fm, which 
was erroneously considered to be small compared to the 
several GeV/fm expected from radiative energy loss. The 
apparent weakness of conventional pQCD collisional en¬ 
ergy loss mechanisms was also supported by parton trans¬ 
port theory results [32|-[3^, which showed that the typ¬ 
ical thermal pQGD elastic cross section, aei ~ 3mb, is 
too small to explain the differential elliptic flow at high 
Pt > 2 GeV and also underestimates the high pT quench¬ 
ing of pions. 

In contrast, Mustafa [2^ found that radiative and elas¬ 
tic average energy losses for heavy quarks were in fact 
comparable over a very wide kinematic range accessible 
at RHIC. In Fig. we confirm Mustafa’s finding and 
extend it to the light quark sector as well. The frac¬ 
tional energy loss, AE/E^ from DGLV radiative for u, c, b 
quarks (solid curves; see also App. IB) is compared to 
TG [131 and BT [2^ estimates of elastic (dashed curves; 
see also App. lA). For light quarks, the elastic energy 
loss decreases more rapidly with energy than radiative 
energy loss, but even at 20 GeV the elastic is only 50% 



FIG. 2: Average AE/E for u,c,b quarks as a func¬ 
tion of E. A Bjorken expanding QGP with path length 
L = 5 fm and initial density fixed by dNg/dy = 1000 
is assumed. The curves are computed with the cou¬ 
pling Qs = 0.3 held fixed. For Debye mass yo oc 
(dNg /the gluon mass is the light quark 

mass is pol^., the charm mass is 1.2 GeV, and the bot¬ 
tom mass is 4.75 GeV. Radiative DGLV first order en¬ 
ergy loss is compared to elastic parton energy loss (in 
TG or BT approximations). The yellow bands provide 
an indication of theoretical uncertainties in the leading 
log approximation to the elastic energy loss. 


smaller than the radiative. 

From Fig. [2] we see that for if > 10 GeV light and 
charm quark jets have elastic energy losses smaller but of 
the same order of magnitude as the inelastic losses. But 
due to the large mass effect [35-[3^,[l3|, both radiative 
and elastic energy losses remain significantly smaller for 
bottom quarks than for light and charm quarks, but the 
elastic energy loss can now be greater than inelastic up 
to ~ 15GeV. We present both TG and BT as a measure 
of the theoretical uncertainties of the Coulomb log (see 
App IC for benchmark numerical examples). These are 
largest for the heaviest b quark. As they are not ultra- 
relativistic, the leading log approximation [ 23 , Hll breaks 
down in the kinematic range accessible at RHIC. More 
rigorous computations of elastic energy loss [1^ and nu¬ 
merical covariant transport techniques [13] can be used to 
reduce the theoretical uncertainties in the elastic energy 
loss effects. 

Theoretical Framework. 

The quenched spectra of partons, hadrons, and leptons 
are calculated as in [Tlj from the generic pQCD convo¬ 
lution 


Ed^cr{e) 

dp^ 


Eid^a{Q) 
dPi 
D{Q 


>P{E, 


Hq) ® f(HQ 


Ef) 

e). 


( 1 ) 


where Q denotes quarks and gluons. For charm and bot¬ 
tom, the initial quark spectrum, Ed^a{Q)/dp^, is com- 
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^ted at next-to-leading order using the code from [2l|, 
1221 : for gluons and light quarks, the initial distributions 
are computed at leading order as in 0. P{E, ^ Ef) is 
the energy loss probability, D{Q ^ Hq) is the fragmen¬ 
tation function of quark Q to hadron EIq, and / {Hq —> e) 
is the decay function of hadron Hq into the observed 
single electron. We use the same mass and factorization 
scales as in and employ the CTEQ5M parton den¬ 
sities [i^ with no intrinsic kx- As in [i^ we neglect 
shadowing of the nuclear parton distribution in this ap¬ 
plication. 

We assume that the final quenched energy, Ef, is large 
enough that the Eikonal approximation can be employed. 
We also assume that in Au-I-Au collisions, the jet frag¬ 
mentation function into hadrons is the same as in e“'"e“ 
collisions. This is expected to be valid in the deconfined 
medium case, where hadronization of Q —> Hq cannot 
occur until the quark emerges from the sQGP. 

The main difference between our previous calcula¬ 
tion m and the present one is the inclusion of two 
new physics components in the energy loss probability 
P{Ei —> Ef). First, P{Ei Ef) is generalized to in¬ 
clude for the first time both elastic and inelastic energy 
loss and their fluctuations. We note that Vitev was 
the first to generalize the GLV formalism to include ini¬ 
tial state elastic energy loss effects in d-|-Au. In this work, 
Eq. (12) extends the formalism to include final state elas¬ 
tic energy loss effects in AA. 

The second major change relative to our previous ap¬ 
plications is that we now take into account geometric 
path length fluctuations as follows: 


P{ Ei El — Arad — Ael) = 


d(j) 

2tt 


(Px± 


■TAA{x±,b) > 


Nbiniby 

Prad {Arad;L{x±,(l>))^ Pel{Aer,L{x±,(j))). (2) 


Here 


L{x±,(j))= J dTpp{x±-\-Tn{4>))/{pp) (3) 

is the locally determined effective path length of the jet 
given its initial production point x± and its initial az¬ 
imuthal direction (f relative to the impact parameter 
plane {x,y). The geometric path averaging used here 
is similar to that used in and by Eskola et al. [d^ . 
However, the inclusion of elastic energy losses together 
with path fluctuations in more realistic geometries was 
not considered before. 

We consider a diffuse Woods-Saxon nuclear density 
profile [d^ . which creates a participant transverse den¬ 
sity, Pp{x±), computed using the Glauber profiles, Ta(x), 
with inelastic cross section a^iN = 42 mb. The bulk 
sQGP transverse density is assumed to be proportional 
to this participant density, and its form is shown (for the 
?/ = 0 slice) in Fig. [3] by the curve labeled pqgp- The 
distribution of initial hard jet production points, af_L, is 



FIG. 3: Transverse coordinate (x, 0) distribution of sur¬ 
viving pt = 15 GeV, Q = g,u,c,b jets moving in 
direction </> = 0 as indicated by the arrows. Units 
are arbitrary for illustration. The transverse (binary 
collision) distribution of initial jet production points, 
pjet(a;, 0), is shown at midrapidity for Au-|-Au collisions 
at b = 2.1 fm. The ratio pq/pjet (see Eq.([8ll) gives the 
local quenching factor including elastic and inelastic en¬ 
ergy loss though the bulk QGP matter distributed as 
PQGp(2:,0). 


assumed on the other hand to be proportional to the bi¬ 
nary collision density, Taa = d"A(x -h bf2)TA(x — h/2). 
This is illustrated in Fig.[3]by the narrower curve labeled 

PJet- _ 

The combination of fluctuating DGLV radiative 
with the new elastic energy losses and fluctuating path 
lengths (via the extra d’^xxdcj) integrations) adds a high 
computational cost to the extended theory specified by 
Eqs. dmi). In this first study with the extended theory, 
we explore the relative order of magnitude of the com¬ 
peting effects by combining two simpler approaches. 

In approach I, we parameterize the heavy quark pQCD 
spectra by a simpler power law, Ed^aQ/d^k oc 
with a slowly varying logarithmic index n = n(pT)- For 
the pure power law case, the partonic modification fac¬ 
tor, Rq = (prior to fragmentation) is 

greatly simplified. This enables us to perform the path 
length fluctuations numerically via 


Rq = 


d<j) 

2tt 


d? 


Xx 


TAA{x±,b) 


Nbin{b) 

J de(I - e)”P(^(e;T(f_L, ?;>)), 


( 4 ) 


where 

^Q(e;T)=y dxPQ^rad{x;L)PQ^eiie - x-,L). (5) 

Both the mean and width of those fractional energy losses 
depend on the local path length. (See App ID for numer- 
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FIG. 4: Distribution of path lengths (given by Eq. (|3|l ) 
traversed by hard scatterers in 0-5% most central 
collisions; the lengths, L{x_l,4'), S'l'e weighted by 
the probability of production and averaged over az¬ 
imuth. An equivalent formulation of Eq. (0 is Rq — 
f dLl/NbindNhin/dL f de(l - L). Since the 

distribution l/NundNun/dL is a purely geometric 
quantity, it is the same for all jet varieties. Also dis¬ 
played are the single, representative pathlengths, Lq, 
used as input in approach II. Note the hierarchy of 
scales with glue requiring the shortest, then charm, light 
quarks, and bottom the longest effective pathlength. 

ical illustrations of Eq. (I5|) for a fixed L = 5 fm light quark 
case.) 

We emphasize, however, that no externally specified a 
priori path length, L, appears in Eq. O; the path lengths 
are allowed to explore the whole geometry. Fig. |4] shows 
the broad distribution of lengths probed by hard partons 
in approach I. 

In the second approach, we determine effective path 
lengths, Lq, for each parton flavor, Q, by varying fixed L 
predictions Rq {pt, L) and comparing to approach I; see 
Fig. [31 In approach II, Rq{pt,L) is calculated directly 
from Eq. (1) with P{Ei Ef) in Eq. (2) replaced by the 
fixed L approximation 

P{Ei Ei — Ael— Arad , L) W 

Prad(^ Arad j L)®Pal{AaV,L). (6) 

Here, jet quenching is performed via two independent 
branching processes in contrast to the additive convo¬ 
lution in Eq. O- For small energy losses the two ap¬ 
proaches are similar. They differ however in the case of 
long path lengths when the probability of complete stop¬ 
ping approaches unity. In the convolution method, the 
probability of e > 1 is interpreted as complete stopping, 
whereas in the branching algorithm the long path length 
case is just highly suppressed. In both cases we take into 
account the small finite probability that the fractional 
energy loss e < 0 due to fluctuations. 


To illustrate the difference in approach II, consider the 
case of power law initial Q distributions as in Eq. ([3]). In 
this case 

R^q^{pt,Lq) = {{l-PQ{LQ)ni-e^Q{LQ)r)AE. 

(7) 

The branching implementation is seen via the product of 
two distinct factors in contrast to the one quenching fac¬ 
tor in Eq. 0. For small both approaches obviously 
give rise to the same Rq = 1 — ri{eQ). 

Due to the high computational cost in approach I, only 
the TG elastic is used for the heavy quarks and only 
BT for light quarks. The Coulomb log uncertainties are 
estimated only in approach II. 

In both approaches, fluctuations of the radiative en¬ 
ergy loss due to gluon number fluctuations are computed 
as discussed in detail in Ref. MM- This involves using 
the DGLV generalization of the GLV opacity expan¬ 
sion Q to heavy quarks. Bjorken longitudinal expan¬ 
sion is taken into account by evaluating the bulk density 
at an average time r = L/2 [Tl|, [s^. For elastic en¬ 
ergy loss, the full fluctuation spectrum is approximated 
here by a Gaussian centered at the average energy loss 
with variance = 2T{AE'^^pt, L)) [s^. In approach 
I the correct, numerically intensive integration through 
the Bjorken expanding medium provides AE^\pt, L). In 
approach II the r = T/2 approximation is again used; 
numerical comparisons show that for L ~ 2 — 7 fm this 
reproduces the full calculation well. Finally, we note that 
we use the additional numerical simplification of keeping 
the strong coupling constant Ug fixed at 0.3. 

Numerical Results: Parton Level 
In Fig. |5l we show the quenching pattern of Q from the 
second approach for a “typical” path length scale L = 5 
fm, similar to that used in previous calculations m- 
The curves show Rq{pt), prior to hadronization, for 
Q = g,u,c,b. The dashed curves show the quenching 
arising from only the DGLV radiative energy loss. The 
solid curves show the full results after including TG elas¬ 
tic as well as DGLV radiative energy loss. Adding elastic 
energy loss is seen to increase the quenching of all flavors 
for fixed path length. Note especially the strong nonlin¬ 
ear increase of the gluon suppression and the factor ~ 2 
increase of the bottom suppression. The curious switch 
of the u and the c quenching reflects the extra valence 
(smaller index n„) contribution to light quarks. 

Fig. O emphasizes the unavoidable result of using a 
fixed, “typical” path length scale, L, in jet tomography: 
the pion and single electron quenching can never be sim¬ 
ilar. If pions were produced only by light quarks and 
electrons only by charm, then we would expect compa¬ 
rable quenching for both. However, contributions from 
highly quenched gluons decrease the pion Raa while 
weakly quenched bottom quarks increase the electron 
Raa- Therefore, in the fixed length scenario, we expect 
a noticeable difference between pion and single electron 
suppression patterns. 
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FIG. 5: Partonic nuclear modification, Rq{pt) via 
Eq.©, for g, u, c, & as a function of pr for fixed L=5 
fm path length and dNg/dy = 1000. Dashed curves 
include only radiative energy loss, while solid curves in¬ 
clude elastic energy loss as well. 

The solid curves of Fig. [3] labeled by the parton fla¬ 
vor Q = g,u,c,b show the relative transverse coordinate 
density of surviving jets defined by 

Pq{x,(I)) = pjet{x) J de(l - e)"'P^(e;L(x, (/>)). (8) 

PQ is given by the initial transverse x production distri¬ 
bution times the quenching factor from that position in 
direction cj) with hnal momentum pp. The case shown is 
for a pt = 15 GeV jet produced initially at (x,0) and 
moving in the direction </> = 0 along the positive x axis. 
The quenching is determined by the participant bulk 
matter along its path pQGp(a; -I- vt, 0, t), and varies with 
X because the local path length L = T(a;,0, 0) changes 
according to Eq. 

What is most striking in Fig. [3] is the hierarchy of Q- 
dependent length scales. No single, representative path 
length can account for the distribution of all flavors. In 
general heavier flavors are less biased toward the surface 
(in direction (j)) than lighter flavors since the energy loss 
decreases with the parton mass. Gluons are more sur¬ 
face biased than light quarks due to their color Gasimir 
enhanced energy loss. In addition, note the surprising 
reversal of the u and c suppressions, also seen in Fig. [S] 
Fig. [3] shows that the energy loss of c is somewhat less 
than for u; however, the higher pp power index, n, of 
c relative to u - as predicted by pQCD and due to the 
valence component of u - compensates by amplifying its 
quenching. 

However, none of the distributions can be categorized 
as surface emission. The characteristic widths of these 
distributions range from Ax « 3 — 6 fm. We show be¬ 
low that such a large dynamic range of path length fluc¬ 
tuations is essential for consistent reproduction of both 


electron and pion data. 

We turn next to Figs. [3]and© which show the inter¬ 
play between the dynamical geometry seen in Fig. ©and 
the elastic-enhanced quenching of partons. In Figs. [3] 
and © the solid green curves labeled “DGLV-I-TG/BT: 
Full Geometry” are the results using approach I based 
on Eq. O- The curves labeled TG and BT are from 
approach II based on Eq. O- The effective hxed Lq in 
II were taken to match approximately the green curves 
in which full path length fluctuations are taken into ac¬ 
count. This procedure is not exact because of the differ¬ 
ent numerical approximations involved, but the trends 
are well reproduced. The Lq are determined only to 
~ 0.5 fm accuracy, as this suffices for our purposes here. 
We show the comparison between approaches I and II for 
heavy quarks in Fig. [3] using = 4.5 and Lf, = 5.5 fm 
and for gluons and light quarks in Fig. ©using Lg = 3.5 
and Lu = 5.0 fm; see Fig. © for a visual comparison of 
the input length distributions used. This hierarchy of Q- 
dependent length scales is in accord with that expected 
from Fig. © 

Note that in comparison to the fixed L = 5 fm case in 
Fig.© geometric fluctuations reduce gluon jet quenching 
in Fig. ©by a factor ~ 2. Nevertheless, even with path 
length fluctuations the gluons are still quenched by a fac¬ 
tor of 10 when elastic energy loss is included in addition 
to radiative. 

The amplified role of elastic energy loss is due to its 



FIG. 6: Heavy quark jet quenching before fragmentation into 
mesons for dNg/dy = 1000. Solid green curves show the re¬ 
sults of approach I based on Eq. © including full geometric 
path length fluctuations and DGLV radiative and TG elastic 
energy loss for c and b quarks. Upper and lower yellow bands 
show predictions using approach II via Eq. 0 with effective 
path lengths taken as Lc = 4.5 and Lb = 5.5 fm. As pre¬ 
viously noted in Fig. [2] the difference between TG and BT 
curves indicates an estimate of the magnitude of theoretical 
uncertainties in the elastic energy loss. 
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FIG. 7; As in Fig.[6]but for light u, d quarks and gluons. 
The yellow bands are computed in this case with effec¬ 
tive g, u path lengths Lg = 3.5 and Lv. = 5.0 fm based 
on Eq. 0. Note that charm and light quark quenching 
are similar in this pr range. 



FIG. 8: The consistency of the extended jet quench¬ 
ing theory is tested by comparing its prediction to 
the nuclear modification of the tt® spectra observed by 
PHENIX [J. 


smaller width for fluctuations relative to radiative fluctu¬ 
ations. Even in moderately opaque media with L/X ~ 10, 
inelastic energy loss fluctuations are large because only a 
few, 2-3, extra gluons are radiated [J. Thus, gluon num¬ 
ber fluctuations, 0{l/y/Ng) lead to a substantial reduc¬ 
tion in the effect of radiative energy loss. On the other 
hand, elastic energy loss fluctuati ons a re controlled by 
collision number fluctuations, 0{^/X/L), which are rela¬ 
tively small in comparison for a significant proportion of 
the length scales probed. Therefore, fluctuations of the 
elastic energy loss do not dilute the suppression of the 
nuclear modification factor as much as Ng fluctuations. 
The increase in the sensitivity of the final quenching level 
to the opacity is a novel and useful byproduct of includ¬ 
ing the elastic channel; see Fig. [Tl]in Appendix D. The 
inclusion of elastic energy loss significantly reduces the 
fragility of pure radiative quenching (45j | and therefore 
increases the sensitivity of jet quenching to the opacity 
of the bulk medium . 

Numerical Results: Pious and Electrons 
We now return to Fig. [T] to discuss the consequence of 
including elastic energy loss of c and b quarks on the 
electron spectrum. The inclusion of the collisional en¬ 
ergy loss significantly improves the comparison between 
theory and the single electron data. That is, the lower 
yellow band can reach below Raa ~ 0.5 in spite of keep¬ 
ing dNg/dy = 1000, consistent with measured multiplic¬ 
ity, and using a conservative as = 0.3. A large source 
of the uncertainty represented by the lower yellow band 
is the modest but poorly determined elastic energy loss, 
AE/E « 0.0—0.1, of bottom quarks (see Fig. [2). There is 
additional uncertainty from the relative contributions to 
electrons from charm and bottom jets. The dashed lines 
show an extreme version of this in which charm jets are 


the only source of electrons. If the charm to bottom ra¬ 
tio given by FONLL calculations is accurate, the current 
data suggests that even the combined radiative-l-elastic 
pQCD mechanism is not sufficient to explain the single 
electron suppression. 

As emphasized in m , any proposed energy loss mech¬ 
anisms must also be checked for consistency with the ex¬ 
tensive pion quenching data [l|, for which preliminary 
data now extend out to pr ~ 20 GeV. This challenge 
is seen clearly in Fig. [U where for fixed L = 5 fm, 
the addition of elastic energy loss would overpredict the 
quenching of pions. However, the simultaneous inclusion 
of path fluctuations leads to a decrease of the mean g 
and u,d path lengths that partially offsets the increased 
energy loss. Therefore, the combined three effects con¬ 
sidered here makes it possible to satisfy Raa < 0-5 i O.I 
without violating the bulk dNg/dy = 1000 entropy con¬ 
straint and without violating the pion quenching con¬ 
straint Kaa ~ 0.2 ±0.1 now observed out to 20 GeV; see 
Fig. [5] We note that the slow rise of Kaa with pp in the 
present calculation is due in part to the neglect of initial 
kp smearing that raises the low pp region and the EMC 
effect that lowers the high pp region (see i). 

Conclusions 

The elastic component of the energy loss cannot be ne¬ 
glected when considering pQCD jet quenching. While 
the results presented in this paper are encouraging, fur¬ 
ther improvements of the jet quenching theory will be 
required before stronger conclusions can be drawn. 

From an experimental perspective, there is at present 
significant disagreement between measured p±p to elec¬ 
tron baselines [3, d] . In addition, direct measurement of 
D spectra will be essential to deconvolute the different 
bottom and charm jet quark dynamics. 
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On the theoretical side, more work is needed to sort out 
coherence and correlation effects between elastic and in¬ 
elastic processes that occur in a finite time and with mul¬ 
tiple collisions. Classical electrodynamics calculations 
presented in suggested that radiative and elastic pro¬ 
cesses could destructively interfere over len gths far longer 
than previously thought. As described in [^, a proper 
accounting of the current shows finite size effects persist 
out only to the expected lengths of order the screening 
scale, 1/^D^ 1 fm- Additionally, work on the quantum 
mechanical treatment of elastic energy loss in a finite 
medium l50,[IH also concluded that finite size effects on 
Raa remain small except in peripheral collisions. 

There are several other open problems that require fur¬ 
ther study. The inclusion of all the initial state effects 
from 0 will be needed to fully check the consistency of 
the pion predictions with the data. Only an approximate 
fluctuation spectrum for elastic energy loss has been in¬ 
cluded here; still needed is an examination of the effect 
of the full fluctuation spectrum. 

The radiative and elastic energy losses depend sensi¬ 
tively on the coupling, oc and oc a^. 

Future calculations will have to relax the current fixed 
as approximation. In [s^ . the running of the coupling 
is seen to increase the magnitude of the elastic energy loss 
and alter the energy dependence. More complete calcu¬ 
lations of both radiative and elastic energy losses will 
involve integrals that probe momentum scales that are 
certainly nonperturbative. Therefore it will be impor¬ 
tant to study the irreducible uncertainty associated with 
the different maximum as cutoff prescriptions commonly 
used. 
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I. APPENDIX 


A. Collisional Energy Loss 

The leading logrithmic expression for the elastic energy 
loss of a jet with color Casimir Cr in an ideal quark-gluon 
plasma wit h rif active quark flavors and temperature, T, 
is given by 

= CR-KalT'^{l + ^)f{v) log(Hc) (9) 
ax D 

where the Coulomb log is controlled by the ratio 
that involves relevant minimum and maximum momen¬ 


tum transfers or impact parameters. For scattering in an 
assumed ultrarelativistic (m = 0) gas of partons the jet 
velocity dependence is 

/(") = ^ (" + - 1) 1 (10) 

Estimates for differ below asymptotic {E ^ T) ener¬ 
gies and are given in [2^, and [1^ that we denote 
by Bj, TG, and BT respectively: 

= {4Ej,T) / 

f (^2^ 0.85 EpT^ / (4) Ep:^^ 

I ^4^ 0.36 / (4) Ep<^^ 

( 11 ) 

with the crossover between Ep ^ ^ and Ep 
being taken at Ep = for numerical computation. 


Hbj 

Htg 

Hbt 


B. DGLV Radiative Energy Loss 


For completeness we also record the DGLV formula for 
radiative energy loss used in our calculations. We neglect 
finite kinematic limits on the momentum transfer q inte¬ 
gral, and perform the finite 0 < fcj_ < kmax = 2px{l — x) 
and azimuthal 0 < (/) < 27r integrations analytically. The 
mean fractional radiative energy loss can be then evalu¬ 
ated numerically from the expression 


AR 


ind 


E 


where 


CpUs L 

TT A„ 


Ep+p 


dx 


4,p-‘q'^dq 


_ Jo (4|E)^ + (g2+^2)2 

x{A\ogB + C) (12) 


/32 = m^{l-x) + M'^ 




^QQ - 


9 ' 

Pg^gg 

97ra^ 


PqO'qg 


and aqg = 


(13) 

(14) 

(15) 


We employ the asymptotic 1-loop transverse gluon 
mass nig = The A,B,C functions denote 

(16) 

h 

^ ^ {(i^ + K){p^Q-g+QtQi + Qtfp) 

/ 32 (/ 32 (q- -K)- QpK + Q+Q+ + fpfp) 

c = 2 q^pf [/3^M^(V-M^) + /3^(/3^-M^)1^ 

+g +(/34 - 2g2Q+) + - 3q^ + 

+ 2 q^Q+) + 3/32 g2Q-] ( 18 ) 
















Radiative 

Gollisional 

AE 

Ejst (GeV) 

DGLV 

Bj 

TG 

BT 

10 

0.2111 

0.2022 

0.1594 

0.1596 

11 

0.2126 

0.1894 

0.1506 

0.1552 

12 

0.2129 

0.1782 

0.1430 

0.1621 

13 

0.2123 

0.1683 

0.1358 

0.1530 

14 

0.2110 

0.1596 

0.1294 

0.1450 

15 

0.2093 

0.1518 

0.1237 

0.1379 


TABLE I: benchmark test cases for a charm jet (m = 1.2 

GeV, Cr = |) with fixed as = 0.3. The path length is L = 5 
fm; he = 0.197 GeV fm. The density is ^ = 1000 (n/ = 0), 
giving T = 0.2403 GeV, /r = 0.4666 GeV^, and \g = 1.2465 
fm. For radiative, the q integration limits are taken to be 
0.0001 and 50. 


where the abbreviated symbols denote 

K = k^max = 2Ta:(l - a;) 

Qt = 

Qk = 

fp = 

with f{x,y,z) = \/x'^ + 2x^y + 


(19) 

( 20 ) 
( 21 ) 
( 22 ) 
(23) 


C. Benchmark Numerical Examples 

In this section we record numerical benchmark cases 
of both the elastic and radiative mean energy loss to il¬ 
lustrate the above formulas. Consider a uniform Bjorken 
cylinder with density 


P{r) 


1 dN 
dy 


(24) 


We assume R = 6 fm. The temperature evolves as 


T(r) 


1.202 9nf + 16) 


where n/ is the number of active quark flavors. The effec¬ 
tive static approximation simulates the effect of Bjorken 
expansion by evaluating T at t = L/2, where L is the jet 
path length to the the cylinder surface. The gluon den¬ 
sity is computed from pg = x 16T^, and the density 
of quarks plus antiquarks is pq = x 9n/T^. The 
Debye mass squared is p? = 47ra^r^(l -|- ^). In Table I 
the results for n/ = 0 are given for a charm jet of energy 
10 < T; < 15. 



FIG. 9: Example of the induced DGLV radiation gluon 
number spectrum for a 15 GeV up quark, with path 
length L = 5 fm, dNg/dy = 1000, Os = 0.3, and 
Ax=118.7 fm^ 


energy loss for a specific case of a 15 GeV up quark jet 
with path length 5 fm. Fig. [9] shows the first order in¬ 
duced gluon number distribution dNg/dx for this case. 
Fig. [TO] shows the fractional radiative energy loss distri¬ 
bution taking into account Poisson fluctuations of the 
gluon number computed as in [^ . The finite probability, 
P{ng = 0) = 0.2377, of radiating zero gluons contributes 
a (5(e) that is not shown. There is also a finite probablity, 
0.0213, of complete stopping with e > 1. 

The width of the elastic energy fluctuations seen in 
Fig. [11] is significantly smaller than the radiative width. 
The narrowing of the convoluted elastic plus radiative 
distributions significantly reduces the distortion effects 
due to fluctuations. Because of the steep px fall off 
of the initial unquenched parton spectra, the smaller 
width of the elastic energy fluctuations considerably am¬ 
plifies the quenching effect due to collisional energy loss 
in comparison to the larger but much broader radiative 
contribution. In terms of an effective renormalization 
(e) ^ Zeffie) as discussed in [^, Zeff is closer to unity 
than the renormalization Zrad ~ 0.5 characteristic of 
pure radiative energy loss distributions. 


D. Energy Loss Fluctuation Spectrum 


This section illustrates the fluctuation spectra of in¬ 
duced gluon number and the distribution of fluctuating 
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FIG. 10: The distribution of fractional radiative en¬ 
ergy loss (or gain), e = AE/E, for the case considered 
in Fig. 1^ The narrow (green) curve corresponds to 
BT elastic energy loss fluctuations; based on the input 
from Fig. [9] the lower, broader (red) curve corresponds 
to inelastic energy loss due to gluon number fluctua¬ 
tions. The part of the radiation spectrum with e > 1 is 
replaced with 0.02135(1 — x), and there is a contribu¬ 
tion from zero gluon number fluctuations, 0.23775(a:), 
not shown. The continuous part of the radiative distri¬ 
bution has integrated norm 0.7374 with mean 0.2026. 
The top (black) curve corresponds to the convolution 
of elastic and inelastic energy fraction fluctuations. 



FIG. 11: The ratio of rms width, o'(L), to the mean 
fraction energy loss (e) for radiative, elastic and con¬ 
volved energy loss distributions is shown as a function 
of the path length, L, for the Bjorken expanding plasma 
with dNg/dy = 1000. The case of an up quark jet with 
E = 15 GeV is shown. Notice that the elastic distri¬ 
bution is significantly narrower than the radiative one. 
This amplifies the effect of elastic energy loss on Raa 
relative to radiative. 
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